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ABSTRACT 


A solution of the Vinti dynamical problem is derived in a 
form that is a clear generalization of the standard form of solu- 
tion of the Kepler problem. Some geometrical results on the 
orbits are given together with a physical model for the potential. 
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THE VINTI DYNAMICAL PROBLEM AND THE GEOPOTENTIAL 

By Diarmuid O'Mathuna 
Electronics Research Center 


SUMMARY 

A new form of the solution of the Vinti dynamical problem 
is derived wherein all three coordinates are expressed in terms 
of one independent variable corresponding to the true anomaly of 
the Kepler problem. 

The relation of the Vinti potential to the geopotential is 
discussed together with the implied inference on the density 
distribution in the geoid. 


INTRODUCTION 

In this report we are primarily concerned with the solution 
of the Vinti dynamical problem, namely, the motion of a particle 
in the field of the Vinti potential. A second concern is the 
implication for geodesy of the close relationship between the 
latter potential and the geopotential. The relevance of the 
dynamical problem to the prediction of orbits of artifical satel- 
lites, first pointed out by Vint’i, is now well recognized. It is 
also of interest that the potential itself suggests a geodetic 
inference on the density distribution in the geoid. 

If the constant c in the Vinti problem is set to zero, it 
reduces to the Kepler problem. Accordingly, motion in the Vinti 
field must be expressible in a form which is a generalization of 
the corresponding form for the Kepler problem. This will be a 
feature of the representation we obtain for the solution, and the 
motivation serves as a guide to our method of obtaining it. 

In the case of the Kepler problem the motion is given in 
terms of the true anomaly (f) by the three formulas for the three 
spherical coordinates, as follows:* 


1 

1 /■, _\ 



- = u = 

— [1 + e cos fj 


(A) 

cos 9 = 

sin i*sin (f + 

w) 

(B) 

tan(cj) - 

fl) = cos i • tan 

(f + (jo) 

(C) 


*An index of symbols appears on pages 82 and 83. 
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in which e denotes the eccentricity, i the angle of inclination, 
and p the semilatus rectum. The angles oj and fi represent the 
argument of perigee and the angle of the ascending node, respec- 
tively. The above system is completed by inclusion of the time- 
angle relationship between time (t) and true anomaly (f) , namely: 


M = n ( t - t ) 


= arctan 


In formula (D) , t Q denotes the time of perigee passage and n, 
the mean motion, is related to the semi-major axis a by the 
relation 



where p is the normalized gravitational constant. The quantities 
p, a, and e are related by the formula 



Our aim here is the derivation of appropriate generalizations of 
the above formulas for the Vinti problem. 

Solutions of the Vinti dynamical problem have already been 
proposed by Vinti and Izsak (refs. 1-6) in which are derived the 
appropriate generalizations of formulas (A) and (B) . A related 
analysis by Aksenov, Gribenikov and Demin has appeared in the 
Russian literature (ref. 7) . In our derivation of the analogues 
of formulas (A) and (B) , we follow a procedure which is substan- 
tially equivalent to that of Izsak; however, the preliminary 
algebraic manipulation is made explicit here and we feel that the 
algebraic relations between the various constants are exhibited 
more clearly. Anticipating the approximation procedure intro- 
duced at a later stage, we note here that from these algebraic 
formulas it is a straightforward matter to make systematic approx- 
imations . 

Corresponding to formula (C) , the third coordinate gives 
rise to an elliptic integral of the third kind and here we have 


e + cos f 


/l - e 2 sin f e Vl - e 2 sin f 


1 + e cos f 
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recourse to an approximation procedure in order to get the de- 
sired representation of the solution. We note that the Vinti 
potential approximates the geopotential only to second order in 
the small parameter? therefore, from the viewpoint of using the 
Vinti model as a basis for Earth-satellite orbit prediction, 
there is nothing lost in making a second-order approximation in 
the solution of the dynamical problem. In deriving the formula 
for the third coordinate such an approximation scheme is adopted 
to give the analog of formula (C) . The same procedure is 
appropriate in deriving the time-angle relationship corresponding 
to formula (D) . It is therefore also permissible to make second- 
order approximations from the algebraic formulas mentioned above 
since doing so is consistent with the overall method. 

It is in dealing with the latter features of the problem — 
namely, in deriving approximate formulas for the third coordinate 
and for the time-angle relationship--that the present treatment 
differs substantially from that of Izsak. In the latter treat- 
ment, besides defining the independent variable f corresponding 
to true anomaly, Iszak further introduces three new auxiliary 
"independent" variables which we here call ^1 , \jj2, and F. The 
last variable corresponds to eccentric anomaly and the first pair 
are the am-functions associated with the Jacobian elliptic func- 
tions which appear in the representations of the first two 
coordinates in terms of f. The third coordinate is then expressed 
in terms of \pi and \p2 ? involving both secular and trigonometric 
terms, and the time-angle relationship is a generalization of 
Kepler's equation which now involves both 4^1 and 4>2 as well as F. 
Both secular and periodic terms appear. The system must then be 
supplemented by an equation giving the relation between 4^1 and 

\p 2 . 


In the present treatment we do not introduce any auxiliary 
variables. We derive the formula relating (p to f in a form that 
is a clear generalization of formula (C) . To complete the 
solution we need the relation between f and t; this is derived 
by similar techniques and appears as a generalization of the 
time-angle relation (D) . Both relations are approximate, valid 
to second order. 

The above considerations occupy the first six sections of 
this report. In Section 1 the problem is formulated. The first 
integrals are derived in Section 2. In Sections 3 and 4 we go 
through the algebraic manipulation leading to the analogue of 
relations (A) and (B) for the first two coordinates. In Section 
5 we derive the formula for the third coordinate corresponding 
to relation (C) . The time-angle relationship corresponding to 
relation (D) is derived in Section 6. 


All the relations involve the Jacobian elliptic functions — 
a fact which helps us write the relations in relatively compact 
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form. However, it is both consistent with the approximations 
already made and expedient as a preparation for numerical cal- 
culations to replace these elliptic functions by their second- 
order approximate truncated trigonometric series representations; 
this is done in Section 7. In Section 8 we include a geometrical 
result which gives some qualitative insight on the orbits. 

As far as we know there has not been proposed a real physi- 
cal situation giving rise to the Vinti potential field. However, 
if we consider how closely the Vinti field can be matched to the 
geopotential, we are led to the construction of a hypothetical 
geoid which would give rise to an external potential field of the 
Vinti type. In Section 9 we take up this question and derive the 
density distribution inside a sphere, which induces an external 
Vinti potential field. More generally, we derive the density 
distribution consistent with an arbitrary rotationally symmetric 
geopotential, the coefficients of which are empirically deter- 
mined. This suggests an acceptable "first-order" hypothesis for 
the density distribution within the Earth. Such an inference on 
the mass distribution suggests a starting point from which further 
refinement may be possible. 


1. FORMULATION OF THE PROBLEM 

We define spheroidal coordinates R, a, and <f> by the rela- 
tions : 


X = 

(k 2 - = 2 ; 

jl/2 

sin a 

COS (f) 

(1.1a) 

y = 

(k 2 + c 2 ; 

j 1/2 

sin a 

sin <p 

(1.1b) 

2 = 

R cos a 




(1.1c) 


from which we note the relations with the spherical coordinates 
r, 0, and <J) , namely: 

r^ = R^ + c^ sin^ a (1.2a) 

r cos 0 = R cos o = z (1.2b) 


The metric coefficients 



for the coordinate system (1.1), 
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are given by 


'll 


2 2 2 
R + c cos a 

2 2 
R + c 


2 2 2 
^22 = R + C COS ° 


(1.3) 


^33 = ( r2 + ° 2 ) Sin2 0 


with 


g ± . = 0 for i ^ j. 


In this coordinate system the Vinti potential has the 
form (ref . 1) : 


^R - c^ c cos aj 

772 — 2 2 \ 

^R + c cos a ) 


in which y is the gravitational constant, c is the constant of 
the spheroidal system, and c^ is arbitrary. Note that when c = 0 
the coordinate system becomes spherical and the potential becomes 
the Kepler potential. When the potential [Eq. (1.4)] is expressed 
in terms of spherical harmonics, the constants c and ci can be so 
chosen that it matches the geopotential up to the second zonal 
harmonic; in fact, by a further modification (ref. 4) ci can be 
adjusted so that Eq. (1.4) matches the geopotential up to the 
third zonal harmonic. The deviation of Eq. (1.4) from the geo- 
potential is then in the higher harmonics and of higher order in 
the small (oblateness) parameter. 

We shall confine our attention to the dominant perturbation 
due to the even zonal harmonics. Since we can then set C]_ = 0, 
the computations are reduced considerably; however, this is not 
necessary for our procedure. Accordingly, we have 


1/ 


yR 

2 2 2 
R + c cos a 


(1.5) 
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The kinetic energy T is given by 


T 




—.2 2 2 
R + c cos 

2 2 
R + c 


]* 2 ♦ l ( 


R + c cos* 


a) o' 


L 1/L2 . 2 \( . 2 \ 1 2 

+ ^-(R + c J^sin a J <J) 

so that the Lagrangian L = T - 1/ is given by 


( 1 . 6 ) 


L = 


[-2 , 2 

2 

R + C 

cos a 

2 

2 

R + 

c 


]* 2 * l ( 


+ ^(r 2 + c 2 )(sin 2 a) 4> 2 


R“ + c“ cos“ a) a* 



2 2 2 * 

R + c cos a 


(1.7) 


When c in Eq. (1.5) is chosen to fit the second-zonal 
harmonic in the geopotential, then the coefficient of the fourth 
zonal harmonic in roughly two-thirds of the corresponding co- 
efficient of the geopotential, and all other harmonics have 
coefficients of smaller order than the corresponding coefficients 
in the geopotential. In fact, the coefficient of P 2 n is of order 
c^ n , whereas none of the known coefficients of the zonal harmonics 
in the geopotential are of smaller order than c^. We shall return 
to this question in Section 9. 


2. SEPARATION: THE FIRST INTEGRALS 

Since <j> does not appear explicitly in the Lagrangian, we 
first utilize the fact that it is an ignorable coordinate. The 
third Lagrangian equation reads : 

^pR 2 + c 2 )(sin 2 = 0 (2.1) 

which immediately yields the integral: 

(r 2 + c 2 )(sin 2 ct) 4> = r (2.2) 
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where A3 is a constant representing the polar component of angu- 
lar momentum. 

We now follow the standard procedure for ignorable coordi- 
nates (ref. 8), that is, we form a new Lagrangian by setting 


k 



where the q r (r = l,...,k) are the ignorable coordinates. We 
find that: 



‘2 2 2 


1 

R + c cos 0 

*2 1/2 . 2 2 \ 

2 

2 2 

R + 2 (R + c cos a) > 


R + c 


+ 

UR 

A 2 

j . 1 


n 2 2 2 

2 /r,2 , 2 \ .2 


R + c cos 0 

z (R + c J sin a 


that is, we have a modified Lagrangian with two degrees of free- 
dom. The modified kinetic and potential energies are 


and 



( 2 . 4 a) 


( 2 . 4 b) 


L = T - V 


( 2 . 5 ) 


To achieve separability, we first write the Lagrangian in 
the form: 
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which, if we set 


R = c sinh £ 
takes the form: 

2 


(2.7) 


L = c 


(sinh 2 5 + cos 2 o)[jj 2 + §o 2 J 

r i 1 

2 2 2 
l_c sinh £ + cos aj 


[( uc sinh 5 + X T Z^l) " T- — h] 


( 2 . 8 ) 


This is now in standard Liouville form. For convenience, in 
subsequent manipulations, we set: 


Q x (£) = c 2 sinh 2 £ , 


r x 3 i l 

V, (£) = - yc sinh £ + — 2 — ' 

1 L 2 cosh* £ J 

Q(£,a) = Q 1 (C) + Q 2 ( 0 ) » 


so that 


■ ♦ ^ 2 ) ' 


v - 


Q 2 (o) = c 2 

2 

cos a 

(2.9a) 

v 2 (a) = f- 

1 

(2.9b) 

2 

cos a 



(2.9c) 

. + V 2> 


(2.10) 
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and 


= q( 

'It 2 

+ ¥ 2 ) 

- V 

(2.11a) 

= Ql 

( 1 t 2 

K2- 

*¥ 2 ) 

- IK + v 2 ) ■ 

(2.11b) 


The Lagrangian equations are: 


d 1 

r9L\ 

9L 

d 1 

r 9L\ 

9L 

dt ' 

^1' 

" 35 ' 

dt' 

k 9cr ' 

3a 


(2 . 12a,b) 


Before integrating these equations, we first derive the 
energy integral. Multiply Eq. (2.12a) by 5 and Eq. (2.12b) by a 
and add to get 


'9? 


W a — i 

fdL) 

) + ° dt 1 

'9a ' 


(2.13) 


which, after rearrangement, gives: 

fi 9L • 9L"| i 3L , • 3L , " 3L , - 3L 

• 5 —r + O — r = 5 + d -s— + t, — r + 0 — r 

L 9 E 9a J H 9a 95 3a 


d 

dt 


dL 

dt 


which integrates to give 


5 + a - L = E (constant) 

9? 3a 


(2.14) 


This is the energy integral. However, when the Lagrangian has 
the form (2.11), we have: 


r OL • 3L 
€ — t + o — 
9? 3a 


s. 9 T . • 3T 

5 — r + a — 7- 

3 5 9 a 


2T 


9 


and the energy integral takes the form 


T + V = E . 


(2.15) 


We now introduce the form (2.11a) into Eqs . (2.12) and get: 





2 ° > 


9V 

3 ? 


(2.16a) 


d 

dt 



dQ 2 

da 



1 • 2 \ 
2 ° ) 


9V 

9a 


(2.16b) 


If we multiply Eq. 


(2.16a) by Q£ , we find 


QK 




+ ¥ 2 





which integrates to give: 


\ Q 2 K 2 = EQ 1 - V x + X x (2.17) 


where is a constant of integration. If we multiply Eq. (2.16b) 
by Qa and proceed in a similar manner, we obtain: 
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(2.18) 


1 ^ 2-2 
2<> a 


eq 2 - v 2 + a 2 


where A. 2 is a second constant of integration. Adding Eqs . (2.17) 

and (2.18), we have: 


QT = Q (E - V) + X 1 + A 2 , (2.19) 


which, in view of the energy integral of Eg. (2.15), implies that 


A 1 + A 2 = 0 * ( 2 . 20 ) 


We shall be primarily concerned with cases with negative energy, 
that is: 


E = -a 2 , (2.21) 


and as a consequence A. 2 must be positive. This fact will be 
evident from Eq. (2.25b) below. Accordingly, we set 


X 


1 




( 2 . 22 ) 


2 2 

and now writing Eqs. (2.17) and (2.18) in terms of a and X , we 
have : 


1 2*2 

JqV 


- T " v ! " a2 °i 


(2.23a) 


U*2 
2° ° 


V " V 2 - 


(2.23b) 


If Eqs. (2.23) are added, we obtain 


T + V = -a , 


(2.24) 
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so that the energy relation is implied by Eqs. (2.23) and so can 
henceforth be ignored. 

We now use relations (2.9) to write Eqs. (2.23) explicitly 
and reintroduce the variable R from (2.7) to get 


•. / r»2 ^ 2 2 \2 „ ,2 A_ 

1 IB— 1_° c °? °.L r 2 = - L. + pE + _1 

2 R + c 2 


2 R 2 + c 2 


2 2 
- a R 


(2.25a) 


l /„2 , 2 2 \2 -2 

^R + c cos aj o = 


2 ■ 2 
z sm o 


2 2 2 
- a c cos o 


(2.25b) 


The pair of Eqs. (2.25), together with Eq. (2.2), constitutes 
the set of equations, the solution of which is the complete solu- 
tion of the dynamical problem. They are the first integrals — 
or action integrals — with constants a 2 , A 3 , and A. The first 
two represent the energy and polar component of angular momentum, 
respectively. The constant A, which also has the dimension of 
angular momentum, does not have an obvious physical interpretation, 
except when c = 0 , in which case it is the magnitude of the angu- 
lar momentum vector. 

We next introduce a true anomaly f defined by 


df 

dt 


A 


2 2 2 
R + c cos 


(2.26) 


where A is a constant, having the dimension of angular momentum, 
which will be defined later; in the Kepler case (c = 0) , the 
angle f has a simple geometrical interpretation. If we let ' 
denote differentiation with respect to f, then Eqs. (2.25) take 
the form 


A 2 R ' 2 = -(r 2 + c 2 ) J \ 2 



sin a 


- 2 jjR + 2a 2 R 2 J + A 2 c 2 

„ 2 2 2 
2a c cos a 


(2.27a) 

(2.27b) 
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and Eq. (2.2) becomes: 


A 3 

2 2 
R + c 

2 

cos o 

A 

(2 2 
(R + c 

• 2 

sm a 

A 3 

1 

c 2 1 

A 

_sin 2 a 

R 2 + c 2 J 


(2.28) 


We now proceed to a consideration of the above equations 
individually . 


3. ADJUSTMENT OF EQUATIONS (2.27) 

In this section, we go through a lengthy procedure of 
algebraic manipulations in order to write Eqs. (2.27) in a form 
which immediately yields the solution we seek. 

The constants appearing in Eq. (2.27) arose naturally in 
the mathematical separation, and are determined readily from 
initial conditions. The constants that will appear in the final 
form are more "natural” in the physical-geometrical sense. To 
determine them in terms of initial conditions, it is necessary 
to determine them in terms of the constants appearing in 
Eq. (2.27). This is the aim of this section. 


The Equation for R 

We first write Eq. (2.27a) in the form: 



We now associate with the energy and angular momentum constants 
two length scales a 0 and p 0 (corresponding to the semi-major axis 
and semi-latus rectum, respectively) by setting 


a 


o 



P« = 



(3.2) 
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and Eq. (3.1) then reads: 


-ii R' 2 = -[(R 2 + c 2 ) (e 2 - 2a 0 R + a 0 p 0 ) - v 2 c 2 a 0 P 0 ] (3.3) 


where the dimensionless parameter v is defined by 


v = 


(3.4) 


We shall find it convenient for the subsequent analysis to 
have Eq. (3.3) written in dimensionless form. To do so we intro- 
duce dimensionless parameters l and n by the relations 


l 


2 

o 




(3.5) 


and define the dimensionless independent variable y by the rela- 
tion 


R 


a 


o 


y • 


Then, in these terms, Eq. (3.3) reads: 


A" 


0 2 2 
2a a 

o 


y 


,2 _ 


= ~[(y 2 + *y)(, 2 -2 y + l 2 o ) - n 2 vV] 
= -\y A - iy 3 + ^(i+nV )*/ 2 

- 2n 2 i*y + n 2 £^(i-v 2 )] 


(3.6) 


Note that in the Kepler case the parameters v and Z Q have a 
geometric interpretation. The ratio of polar angular momentum 
to "total" angular momentum is measured by v which corresponds, 
therefore, to cos i where i is the inclination in the Kepler case. 
The parameter £ Q measures the ratio of " latus rectum" to "major 
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2 

axis" and so corresponds to (1-e^) where e^ is the Kepler 
eccentricity. The small parameter n measures the "oblateness" 

(or, alternatively, the inhomogeneity, c.f. Section 9) against 
the characteristic dimension of the orbit. We shall ultimately 
approximate the solution in terms of this small parameter. We 
first aim at a form of the solution suitable for making systematic 
approximations . 


Once the initial conditions have been specified, the 
quantities A 3 , A and are known from Eqs . (2.2) and (2.5) and 

so the length scales a Q and p Q are immediately determined. We 
therefore refer to a Q and p 0 as the fundamental length scales. 
Similarly, we shall refer to v, £ 0 , and n as the fundamental 
parameters. Most of the subsequent algebraic manipulation is aimed 
at expressing the constants which shall appear in the representa- 
tion of the solution in terms of these fundamental quantities. 


We start by decomposing the quartic on the right side of 
Eq. (3.6) into two quadratic factors such that each factor is a 
"perturbation" of the corresponding factor in the Kepler case: 
we set 


2 ,, 2 . 2 , 
w - 2ri £ h u 
J o o 


l 


• hr - 2 


yy 

- 2 r 1 2 ^£/ + n 2 £g(l - V 

^ [ 

s - h (l - ri^&^h \ 1 
0 o\ 0 0/ J 

s 

0 

2 2 

1 - n a h 

00 J 

h o)^ + 

^-n 2 *X) 2 ] • 


(3.7) 


With this choice the identity requirements on the first, second, 
and fourth coefficients are automatically satisfied. The identi- 
fication of the third and fifth coefficients gives the pair of 
equations for the determination of h and s , namely: 



2.2 

n 




t \ 2 1 


2 \ 2 

h ) 

o 0/ J 


2 2 

1 - ^*o h o 


+ 4n^h (l - t) 2 Z 2 h ) = 1 + r\ 2 1 2 
o \ 00/ o 


(3.8a) 
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and 


U - " 2<! o h o)[ s ° - “of 1 - " 2 *o h o) 2 ] 

2 

s o 

which after some rearrangement may be written: 



(3.8b) 


(3.9a) 


(3.9b) 


This pair of algebraic equations can be uncoupled if we associate 
with s Q the related quantity q Q defined by 


s 


o 



nVh V 

O O' 


(3.10) 


Inserting Eq. (3.10) into Eq. (3.9b) and rearranging, we get h Q 
in terms of q 0 in the form 



(3.11a) 


and hence 


1 - 


2 2 

n rh 


l - 


2 . 2 
n H q 
o M o 


o o 


|\ 2 . 2 / , ‘ 2 \ 21 

} - n ^ol 1 " V K 


(3.11b) 
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Substituting s Q and h Q from Eqs. (3.10) and (3.11) into 
Eq. (3.9a), we get the equation for q Q , namely: 


[l + nV( 1 - v 2 )q 2 ][l - n 2 l 2 (l - v 2 )q 2 ] 2 
+ 4n 2 q 2 (i - n 2 n 2 q 0 )[l - q D (l - v 2 )] 

- iJ 1 * i 2 *o)[l - r ' 2 e o ( 1 - v2 )%] 2 • (3 ' 12) 


This sixth-order equation cannot be solved algebraically. How- 
ever, if we note that the root we seek is that which tends to 1 
as n tends to zero so that 


q o = 1 + o(n 2 ) , 


(3.13) 


then we see that the determination of the root as a power series 
in r|2 i s a straightforward matter. In fact, if we anticipate the 
later approximations when effects of order ri b are neglected, it 
would be consistent to replace Eq. (3.12) by the cubic equation: 


[l - >A 4 (l - v 2 ) 2 ][l - n 2 * 2 (l - v 2 )q 2 ] 

+ 4 ’1 2 ( 1 - ^ohot 1 - lot 1 - v2 )] 

= + " 2i!2 )[ 1 + - v2 ) - - V 2 ), 2 ] 

(3.12*) 

which can be solved algebraically. However, it is probably more 
economical to determine the root as a power series in ri^. 

With q Q determined from Eq. (3.12), h Q and s^ follow 
immediately from Eqs. (3.11a) and (3.10), respectively, and so 
the coefficients in the quadratic on the right of Eq. (3.7) are 
all determined. From now on we shall consider s Q , q Q , and h Q 
known in terms of the fundamental parameters. We add the remark 
that taking Eqs. (3.13) with (3.10) and (3.11) shows that 
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Iq = 1 + o (n 2 ) , h o = V 2 [i + o(n 2 )] . 


(3.14) 


We now combine Eqs . (3.6) and (3.7) and factor out the con- 

stant term in the second factor. If we note from Eq. (3.10) that 


2 2 

l - n £ h 

o o 


= q (l - r| 2 & 2 h ) = -»/s q 

^o\ o o/ V o^o 


(3.15) 


and substitute for a from Eq. (3.2), we get: 


* 2 ,2 

A y = 


[ 


ua £ 2 T -t/s q 

K o o L „ V o^o , 
11 - 2 3 y + 

jr 


^O 2 ] 

J 


2 „ 2 n 2,, , 2 „ 4 ( q o h o) 

V -2r|£h y + r\ l _ ■ ■ ■ — 

°° ° \/q s 

" M o o 


(3.16) 


where we have interchanged the order of the factors. Written in 
terms of the original variable R and again noting Eqs. (3.2) and 
(3.5), Eq. (3.16) becomes 


2 2 
A^R' ^ = 


r 




R + 


i^R 2 l 

1 P 

0^0 J 


r 2 - 2n 2 h p r + n 

0^0 


2 ( q o h o) 2 

r — p o 

V q o s o 


(3.17) 


At this point, still following the procedure for the Kepler 
problem, we introduce a new independent variable u, defined by 


u 


1 

R 


so that 


u' 



(3.18) 
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and Eq. (3.17) takes the form: 


A 


4 „ 2 - 2 S u + JoJ 

q o L Po a o p oJ 


i 0 2. , 2 ( q o h o) 2 2 

1 - 2ri h Q p o u + n _ — p^u 


Vv 




, 0 2. 2 q o h o) 2 2 

. - Vo“ + q V j 


(3.19) 


This suggests defining the quantity e (corresponding to the 
eccentricity) by the relation ° 


o A/ 

e 2 = 1 - — 
o s 

o 


(3.20) 


and Eq. (3.19) then reads: 


A 2 u ' 2 = hi 

<3^ 


Vv 


2 ~\ 


q s_ — 7T - u - 
° ° \ P. 


. 0 2. , 2 q o “ h o) 2 : 

l-2rihpu + n — p u 

o o 


Vv^ 


The next step is to set 


(3.21) 


Wo (, \ 

u = — (l + e Q w) 


(3.22) 
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which, when introduced into Eq. (3.21), gives as the equation 
for w: 



20 



so that the above equation for w reads: 


,2 ,2 ,2 
Aw' = A 


[ X - " 2 %*o( 3h o - <0] j- x . w 2j 



2 2 2 2 2 
2n + n g^e^w 



(3.25) 


To set this equation up for a straightforward solution, we 
make a final adjustment to the quadratics on the right of 
Eq. (3.25). We write: 


1 - w 2 = J 2 [(l - 6w) 2 - (w - 6) 2 J 
= J 2 (l - 6 2 ) (l - w 2 ) , 


(3.26) 


so that 


J 2 (l - 6 2 ) = 1 . 


(3.27) 


Seeking a similar decomposition for the second quadratic, we 
write : 


1 


- 2 . 2 2 2 2 

2n h-L^QW + n 9 1 e o w 




+ B (w 



(3.28) 


which yields the following relations 
J 2 (a + B<5 2 ) = 1 
J 2 (b + AS 2 ) = n 2 e 2 g 2 
J 2 6 (a + b) = n 2 e Q h 1 


"N 

► (3.29) 


from which, together, with Eq. (3.27), we derive the following: 


i 
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A 


B = 1 


(3.30a) 


2 2 2 
“ n e 0 g l 


(l + 6 2 ) 

(l " <$ 2 ) 


(A + b) 


l + n 


2 2 2 

e o g i 


1 - 6 ‘ 


(a + b) = 


2 u 

n e Q h 1 


(3.30b) 

(3.30c) 


Combining Eqs . (3.30a) and (3.30b), we have the equation for 

6, namely: 


2 u 

n e Q h 1 


. . *2 . . 2 2 2 ' 
1+6 1 + n e Q g^ 


(3.31) 


so that 6 = 0 (n 2 ) • We set 


6 = n e d , 
o o ' 


h = 


1 + 


2 2 2 
n e Qgi 


(3.32) 


and then Eq. (3.31) is equivalent to the quadratic equation for 
d 0 : 

n 4 e 2 hd 2 - d + h = 0 (3.33) 

o o o 


with solution 


d = 
o 


T, ~ 4 2. 2 
- y 1 - 4r) e Q h 

2n e 0 h 


" h [ 


i 4 2. 2 . ./ 4 2\ 2 *I 
1 + n e Q h + 0(n e o ) J 


(3.34) 
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If we substitute for 6 from Eq. (3.32) and use Eq. (3.31) in 
Eq. (3.30b), we obtain: 

A -f B = 1 + n 2 e 2 g l - 2n 4 e 2 hd 0 (3 . 35) 

~ 1 + ^e 2 (g 2 - 2n 2 h 2 ) . (3.35*) 

If we now combine Eq. (3.30a) with Eq. (3.35), we obtain for A 
and B: 


A = 1 - n 4 e Q hd o , B = n 2 e 2 (g 2 - 0 2 hd Q ) 


(3.36) 


Returning to Eq. (3.25) and using Eqs . (3.26) and (3.28), we 

have : 


a 2 w ,2 . z f 1 ~ - g p ) 

• J 4 [(1 - 6w) 2 - (w - 6) 2 ] 

• |a( 1 - 6w) 2 + B (w - 6 ) 2 J , (3.37) 

4 4 

or dividing across by J (1 - 6w) , we have 


r w_. i\ A 1 

- n2 v^( 3h o - O] 

Lj 2 (1 - 6w)J 

q o 


• [> • (f^fl 
■ [» * 
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We therefore write: 



w - 6 v + 6 

v ~ 1 - 6w or W 1 + 6v 

/ 

(3 

so that. 

again using Eq. (3.27), we have: 




i w' 

v 2 ' 2 
3 (1 - 6w) z 


(3 

and Eq . 

(3.38) takes the form: 

o r o . / > n 



> 

tvj 

< 

.2 X I 1 ‘ " VVo( 3h o - %)}L 

% L 

- V 2 ][a + BV 2 ] 



, - ^Vv^K - %)] ( , + B j 

q o 


introducing (A + B) and B from Eqs . (3.25) and (3.36), respec 

tively, into Eq. (3.40) gives: 

a2v ' 2 = ^I 1 ‘ n2 Vv^( 3h ° ‘ q °)] 

•[l + r , 2 e 2 (g 2 . 2n 2 hd Q )][l - v 2 ] 

[* + " 2 <s o( g i ■ 2 r| 2 hd o)] 


39a) 


39b) 


.40) 


.41) 
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where, by Eqs . (3.24), (3.32), and (3. 34), the constants h, 

and d Q are determined in terms of the initial parameters. 

Next we do a corresponding analysis on the equation for a. 


The Equation for a 
Rewriting Eq. (2.27b), we have 


.2 ,2 

A a ' 


. 2 

sin a 


2 2 „ „ 

3 _a c . 2 2 

— 77 ~ 2 — sin a cos a 

A^ A^ 


. 2 

sin a 


] 

] 

[(l-vV (l+ ^)oo.*o + cos 4 aj 


r . 2 

sm a - 

[ si 


A“ I . 2 2 .22.2 2 

— 2 — sm a - v - £ o n sm a cos a 

sin z a " 


(3.42) 


where we have substituted from Eqs. (3.2), (3.4), and (3.5). 

We now define a parameter m D in terms of the fundamental 
parameters 1 , v, and n by the relation: 


2m 

o 


1 + 


- - 2 ) 
(1 + ^ 2 ) 2 


1/2 


so that 


m Q = 1 + 0 ( n 2 ) . 


(3.43) 


(3.43*) 


If we consider the quantity on the right side of Eq. (3.42) as a 
quadratic in cos^ a, its roots are neatly expressed in terms of 
m 0 . In fact, it can be readily checked that we can write 
Eq. (3.42) in the form: 
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.2 ,2 A 

A a' = 


• 2 „ 

sm a 


(l + £ 2 r| 2 )m - n 2 &^cos 2 al 

V 0/0 O J 


1 - v 2 

- cos a 


(i - # 2 ) 


m 


X 2 (l + n 2 ^)m 0 


. 2 

sin a 


1 - 


2 „ 2 
n JL q 


(l * n 2 #» 0 


2 

cos a 


(1 i >n 2 K 


- cos a 


(3.44) 


We now make the substitution: 


cos 



so that 


(3.45) 


-sin O'O' 



(3.46) 


and Eq. (3.44) takes the form: 


*V 


2 





" 2 **(i-» 2 ) .2 

(1 + „ 2 , 2 )V ; 


(3.47) 


where we have reversed the order of the factors. Equation (3.47) 
is what suggests our choice for A. 
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Definition of A 


1 


The only requirement on the quantity A is that it have the 
dimension of angular momentum. As suggested by Eq. (3.47), we 
define A by making the identification: 


K 2 = * 2 (l + nV)m o 


and the pair of equations for v and z, [(3.41) and (3.47)] take 
the form: 


\} n2 Vq 0 5 o( 3h o ~ + n2e o(gi ~ 2n2hd o)] 

q m (l + r] 2 Z 2 ) 

^o o\ o/ 


[> - 


am 

t * 2 2/ 2 0 2. . \ V ) 

1 + n e () (g 1 - 2r\ hd Q ) 


,2 _ 


t - ‘1 


, ^ .2 
(1 + 


The next step is to obtain explicit representations for the 
solutions of the above equations . 


4. EXPLICIT SOLUTION FOR R AND a 


Considering Eqs . (3.49) and (3.50), we set 


[i - n 2 - <a 0 )][i + i 2e o(9i - 2 i 2 ha 0 )] 

Vol 1 + 1 2 t o ) 


» 2 4(* 2 1 - 2 * 2 ™ o ) 

1 + n 2 e 2 (g 2 - 2n 2 hcl 0 ) 


n 2 ^(l - v 2 ) 

( 2 ♦ A 2 ) 2 , 2 


and the pair of equations takes the form: 


V 2 = j 2 [l - v 2 ][l - k 2 (l - v 2 )] 

?’ 2 = [l - £ 2 ](l - k 2 ?2 ] 

with solutions: 

v = cn [^l( f + w l) ,k l] 

Z, = sn ff + w 2 '^ 2 ! 


(4.3a) 

(4.3b) 


(4.4a) 

(4.4b) 


respectively, where and w 2 are arbitrary constants introduced 
by the integration. 

If we define "perigee" as the points at which v' = 0 and 
v" > 0, and if we make the "angle" f have origin at perigee, then 
it follows immediately from Eqs. (4.3a) and (4.4a) that 


Wl = 0 . (4.5) 

If we define the "angle of perigee" w as minus the value of f at 
the first equatorial crossing (? = 0) , then from Eq. (4.4b) we 
have 


and the solution [Eq. (4.4)] is now 

v = cn|j^f,k^l , Z = sn |f + to,k 2 J • (4.7a) 


It remains to express R and a in terms of f. 

For the R expression we first take Eq. (4.7a) and insert it 
into Eq. (3.39a). Then the resulting expression for w is sub- 
stituted into Eq. (3.22) to give: 
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= U = 


_Vi 


o s o 


1 + e 


cn [^l f,k l] + 6 
° 1 + 6 cn[j 1 f,k 1 ] 


. V^l 1 + 3e 0 ) 


1 + 


/ e 0 + 6 \ 

\1 + 6e Q / cn [^ 1 f,k 1 ] 


1 + 5 cn 


Pi £ ' k i] 


(4.8) 


In a similar manner, if we insert Eq. (4.7b) into Eq. (3.45), we 
have 


cos a = 


/ 1 

2 

- V 

f(l + 



(4.9) 


We now see how to define the final constants of the problem, 
namely those appearing in the solutions [Eqs. (4.8) and (4.9)]. 

We shall refer to these constants as the semilatus rectum p, the 
eccentricity e, and the inclination parameter N, defined by 


P = 


+ 6e o) V q o S o i 1 + n2e o d o) 


(4.10a) 


e + 
« - o 

4 . ^ol 1 

+ " 2<3 o) 

1 + 5e o (l + 

1 - H 2 = 1 V* 2 


( 

. 2„2> 
l + n & , 

^ OJ 

1 o 

solutions [Eqs. 

(4.8) and (4.9)] 

1 1 
R U P 

1 + e cn 

[j 1 £ ,k l] 

1 4- S cn 



(4.10b) 


(4.10c) 


(4.11a) 
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cos a 


Vl - N 2 sn [f + w,k 2 J 


(4.11b) 


We shall find useful the following alternative expressions 
for k^, k 2 , and 6, namely: 


k £i <-> *-> 

, = n e • 


2 2 ( 1+n e O d o) g l “ 2n hd 0 


(i + n 2 d 0 ) 2 1 + n e o( g i" 2n hd o) 


2 2 2 
= n e g 


(4.12a) 


*2 - ' » 2 ) 


(l + n 2 ^)m 0 


= n 2 (i - n 2 ) 5- 2 


(4.12b) 




6 = rTe-d 9 W — rr = 0 e.d 

l + n 2 d o 


(4.12c) 


where the quantities g, £, and d are defined by the above rela- 
tions . 

Note the phenomenon of perigee precession is here indicated 
by the fact that the right-hand side of Eq. (4.11b) has period 
4K 2 7 ^ 2 tt where 


2K 


2 


7T 




(4.13) 


The deviation of the right-hand side from unity measures the 
perigee precession rate. 


5. THE INTEGRATION OF THE <j> EQUATION 

Turning to Eq. (2.28), we substitute for A from Eq. (3.48) 
and noting Eq. (3.4), we obtain (setting R = 1/u) : 
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(5.1) 


<j)' = 


[C 1 + ^oK] 


1/2 


r _^ am 

.2 n . 2 2 * 
Ism a 1 + c u J 


Using relations (3.5) and (4.10), we rewrite the above equation 
in the form: 


<J>' = _ - 


072 


[v 2 + (l + nV)m 0 - l] 

r , 2 q o s o ( 1 + r | 2 e o d o ) 2 ( pu ) 2 


N 

. 2 

sin a 


- N' 


1 + ^Vol 1 + n2e o d o) 2 (P u ) 


(5.2) 


It is now convenient to introduce the symbols: 


= 


_ [v 2 + (l + nV)m 0 - l] 


1/2 


= 1 + 


V 


°(n 2 ) 


(5.3a) 


S 1 = 


= q oM ♦ = 1 + o(n 2 ) 


(5.3b) 


so that Eq. (5.2) reads: 


N 


q l*' = 2 


- N ' 


sm a 


2 , x 2 
n s 1 (pu ) 

1 + r) 2 s 1 (pu) 2 


(5.4) 


or, anticipating the approximation procedure, we write 


= 3 Ns 1 |"r| 2 (pu) 2 - n 4 (pu) 4 l + o(n 6 ) . (5.5) 

sin o L -* 
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In the subsequent computation, the evaluation of each term must 
be done separately. We therefore introduce the auxiliary 
quantities <j>., <J> 2 , and <j> 2 as follows: 



1 


N 


2 

cos 


a 


(5.6a) 


<J>2 = (P u ) 2 / 



(5 . 6b, c) 


so that Eq. (5.5) reads: 


q j_<P ' 


< P ± ~ Ns l( r ' 2< ^ ) 2 “ ^^ 3 ) + °( n6 ) 


(5.7) 


which integrates as 


q-^ - = 4>j_ " Ns i f i 2 ( < t ) 2 ” n 2 <l> 3 ) + o(n 6 ) 


(5.7*) 


where is the constant introduced by the integration. 

We now consider Eqs . (5.6) individually. Starting with 

Eq. (5.6a), we introduce cos a from Eq. (4.11b). If we write: 


F 2 = f + w , 


(5.8) 


we have 


d( ^l _ N__ 

dF 2 “ 1 - (l - N 2 )sn 2 (F 2 ,k 2 ) 

= N nc 2 (F 2 ,k 2 ) 

1 + N 2 sc 2 ^F 2 ,k 2 ) 


(5.9) 
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and so 


dcj> 1 N nc 2 (F 2 ,k 2 )dn(F 2 ^ 2 ) 

^ + » 2 '.c 2 (F 2 . k2 ) 

+ N n ° 2 ( F 2/ k 2 ) t 1 - dn ( F 2’ k 2 )] 

1 + N 2 sc 2 F 2 ,k 2 


N nc 2 1 

( F 2' k 2) 

|dn| 

( F 2' k 2 ) 

1 + 

. t2 2 , 
N sc | 

( F ; 

2 ,k 2 ) 

1 


+ k 2 N sc2 ( F 2- k 2) 

2 [l + N 2 so 2 (F 2 ,k 2 )][l + an(F 2 ,k 2 )] 

N nc 2 (F 2 ,k 2 )dn(F 2 ,k 2 ) 

1 + N 2 so 2 (F 2 ,k 2 ) 


2 , 
nc | 

( F 2' k 2 ) 

1 

_ 1 

1 

2 

1 + N 

2 , 
sc | 

( F ; 

>' k 2 ) 

1 

1 + dn| 

( F 2' k 2 ) 


N nc 2 | 

! F 2' k 2 ) 

|dn| 

( F 2 ,k 2 ) 

1 + 

N 2 sc 2 (F 2 ,k 2 ] 

1 


2 , 2 K _ M 1 1 

L dF 2 J 1 + dn(F 2/ k 2 ) 


where we have substituted from relation (4.1 2 b) and utilized 
Eq. (5.9). 
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For the integration of Eq. (5.10), we first make some 
observations. In connection with the second term on the right, 
we note that 


1 + dn(F 2 ,k 2 ) 


2 + 1 ~ dn ( F 2' k 2 ) 

1 + dn(F 2 ,k 2 ) 


(5.11a) 


1 + 


sn 2 (F 2 ,k 2 ) 


[ 


1 + dn( 


‘( F 2' k 2 )]' 

while for the first term on the right we note that 

N nc 2 (F 2 ,k 2 ) dn(F 2 ,k 2 ) 

( F 2 ,k 2) 


(5.11b) 


, r i "I N nc 

gp-jtan" [Nso(F 2 ,k 2 )]j = — - 


m 2 2 
N sc 


(5.12) 


so that Eq. (5.10) may be written in the form: 

afjjh sc ( F 2' k 2 )] 

.2,2 /a*, 


_ ri Z 

2 VdF, 


- N 


i + k: 


( F 2 ,k 2 ) 


[l + dn(F 2 ,k 2 )] 2 _ 


(5.13) 


Again, with a view to integrating the right-hand side, we note: 

2 2/ s 2 N sc2 

k 2 sn ( F 2' k 2)’ — = k 


dF. 


N sc^(F 2 ,k 2 ) 

2 1 + N 2 sc 2 (F 2 ,k 2 ) 


k 2 N 


2 

nc 

i F 2 ,k 2 ) 

-1 

"2 


2 


_L 

1 - N 

1 

+ N 

sc ( F 2' k 2 ) 



= n MdFT 


- N (5.14) 
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where we have again used Eq. (5.9) and the trick used in arriving 
at Eq. (5.10). Inserting Eq. (5.14) into Eq. (5.13), we obtain: 



The procedure separates out the dominant contribution and it 
appears that it can be continued indefinitely. However, in our 
approximation, we do not seek accuracy beyond the second power of 
an d so we terminate the manipulation at this point. We note 
that* 

dn ( F 2' k 2) = 1 + °( n2 ) 

sn(F 2 ,k 2 ) = sin G 2 + 0 (n 2 ) 

where 

°2 “ 2% F 2 = F 2 [ 1 - T - TT + °( k 2 )] 

If we now carry the approximation to the second power in r\^ , it 
is consistent to replace Eq. (5.15) by 


(5.16a) 

(5.16b) 

(5.16c) 



if-) - tan_1 [ N sc ( F 2' k 2 )] 

- ' 3^N(l-N 2 )sin 2 G 2 
“ - 2^n(i - N 2 ) (l - COS 2G 2 ) 


(5.17) 


*Approximations such as these used in sections 5 and 6 are derived 
in section 7 . 
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or , on rearrangement : 



which integrates to give: 


(5.18) 


— !L ¥~) = tan 1 [ N so( F 2 ,k 2 )] 

- nV N |; +n V(3. N 2 ) ] F2 

4 4 

+ ^|-n(i - N 2 ) sin 2G 2 . (5.19) 


In Eq. (5.19) we have not included the arbitrary additive constant 
since we can consider it already absorbed in the constant 
appearing in Eq. (3.7*). 

We next consider Eq. (5.6b). If we introduce u from 
Eq. (4.11a) and set 


F i - h f 


r 


we have : 


d^2 

1 

+ 

e 

cn ( F i 

- k i) 

dF l 

r 

+ 

6 

cn ( F i 

- k i)_ 


(5.20) 


(5.21) 


Since in Eq. (5.7*) the factor <J > 2 is multiplied by n 2 , it is 
sufficient for our approximation to compute <j> 2 up to the first 
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power in ri . Accordingly, using relation (4.12c), we note that 


1 + e cn(F^,kjJ 
1 + 6 cn^F^,k.jJ 


= 1 + e (l - r| 2 d) cn^F^ ,k^) 


- r| 2 de 2 cn^F^,kjJ + o(n^) 


and so, inserting in Eq. (5.21), we have: 
d<J) 


j-L r- = 1 + 2e (l - nd) cn^,^) + e 2 (l - 4r]d) cn^F^k^ 
- 2r| 2 de 3 cn 3 (F^,kjJ + o(r| 4 ) 


= 1 + 2e(l - n 2 d) 


k 2 

S " 2 ( F 1 

' k ll 

L 

) cn 

( F i ,k i) 

k l 

1 + 

dn| 

( F 1 

- k i) 

cn | 

( 2F l' k l 

) + 

dn 

( 2F i' k ii 


cn ( F l' k l) dn ( F l' k l) 

+ e 2 (l - 4n 2 d) 


1 + dn(2F^,k^J 


- 2n 2 de 3 cn 3 (F 1 ,k 1 ) (5.22) 


4 

where we have now omitted terms of order n . The first bracket 
in Eq. (5.22) was obtained by a manipulation similar to that done 
below Eq. (5.9), while the second bracket was obtained by 
expressing cn^ in terms of the elliptic functions of the double 
argument . 
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In the first bracket in Eq. (5.22) we set 


1 + dn^F 1 ,k 1 ^ 


- I ♦ °(* 2 ) . 


(5.23a) 


and in the second bracket we set: 


1 + dn 


2 

[ ■ 3k^ 

dn ( 2p 1 ' k 1 ) + -4 i =n 2 (2 Fl ,ki) + 0(^ 4 ) 

. _ _ 


(5.23b) 


and we can replace Eq. (5.22) by 


d<j>, 

i i 

3 1 dFt 


1 + 2 e(l - n 2 d) ^cn(F 1 ,k 1 )dn(F 1 / k 1 ) 

2 

^1 2 ( ^ t \ /V, , \\ , e(l - 4 n^d) 

- T sn ( F l' k l) cn ( F l' k l)J + ~2 

• Tl + cn ^ 2 F^ ,k^) dn ( 2 F^ ,k^J - k 2 sn 2 ( 2 F^,k^) 

•[ F - |( on ( 2 F l- k l) + dn ( 2 F l' k l))]] 

- 2 n 2 de 2 cn d (p^,k^) 

1 + f (l - 4 n 2 d) + ^[ 2.(1 - n 2 d)sn(F 1 ,k 1 ) 
+ |-(l - 4 n 2 d) sn ( 2 F ;i ,k 1 )J - k 2 e 

* [( X " n 2 d)sn 2 (F 1 ,k 1 )cn(F 1 ,k 1 ) 

+ j(l - 4 n 2 d) sn 2 ( 2 F 1 ,k 1 ) [l - |-(cn ( 2 F 1 ,k 1 ) 

+ dn^F^k^jJ - 2 r) 2 de 3 cn^F^k-J 


(5.24a) 


(5.24b) 


38 



In the latter terms of the above expansion (namely, the terms 
with kj- as a factor) we now take: 


dn(F 1 ,k 1 ) = dn(2F 1 ,k 1 ) = 1 + o(n 2 ) 


(5.25a) 


sn(F 1 ,k 1 ^ = sin G-j^ + o(r| 2 ) , 
sn^2F^,k^^ = sin 2G^ + o(r| 2 ) 


/ , \ 

„ / 2\ 

cn ( F l' k l) 

= cos G^ + 0 ^ri J 


cn^2F^,k^^ = cos 2G^ + o(n 2 ) 

where 


7T 

r* — 11 Tn — rr» 

r k 2 

i _ 21 _ 

5k 2 

+ o /^k ^ H 

G 1 2K 1 F 1 F 1 

4 

64 

+ 0^k 1 M 


| (5.25b) 


| (5.25c) 


(5 . 25d) 


and also, noting relations (4.12), it is consistent with the 
approximation to replace Eq. (5.24) by 


h dff " 1 + K 1 ' 4r,2d ) + df^ 2 ^ 1 ' •n 2 d)sn(F 1 ,k 1 ) 
+ |(l - 4n 2 d)sn (2F 1 ,k 1 )J - n 2 e 3 



3 

8 


. 2 
sm 


2G^ cos 2G^ 


+ 



cos 




cos 


G 1 + cos 
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and 


d<j>, 

• 4 

D 1 df\ 


r o / 2 2 

l + |- 2n 2 s ( d -2_|_ 


)J +a ^[ 2e ( 1_ r > 2 ‘ a ) sn ( F l' k l) 

|(l-4n 2 d)sn(2F 1 ,k 1 )]- 


2 3 
n e 


3 3 

^ - 2 /sin G 1 sin 2G^ sin 4G^ 

af-19 


16 


64 


,/sin 3G, \ 

I 3 ~ + 3sin G 1 < 5 


We next note from Eqs. (4.12), (3.34), (3.32) and (3.24) that 


d = 2n 2 - l + o(n 2 ) , g 2 = l - n 2 + o(n 2 ) , 

so that integrating Eq. (5.26) and rearranging, we get: 


(5 


1 

r 

2 ^ / 

2 \“|" 

11 

CM 

-e- 

1 + | + 2ne ] 

L + h - 2N ( ; 

L + h)\ 


+ 2e 


[l + n 2 (l - 2N 2 )_ 


( F l ,k l) 

~ X 

+ e [l 4- 4n 2 (1 r 2N 2 )] s „(2F 1 ,k 1 ) 

[ 


2 3 | 7 - 13N 2 „ 3(1 - N 2 ) 

- n e | ^ sin G x - 2 sin 2G 1 

2 2 2 
+ — - j 2 - N ~ sin 3G ± - - 1 - -^ 4 N sin 4G 1 + 1 sin 6G 


■] 


(5 


.26) 


. 27) 


.28) 
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where we have used the additive constant to replace f by F 2 in 
the secular term. The reason for doing this will appear later. 

It remains to calculate <£ 3 . If we substitute for u from 
Eq. (4.11a) into Eq. (5.6c) we have: 


*3 


1 + e cn^f,]^) 
1+6 cn^f,]^) 


4 


(5.29) 


We note that in Eq. (5.7*) cj> 3 has a factor g 4 . Accordingly, in 
the calculation of $ 3 , we may neglect terms of order r)2 # For Q ur 
approximation, therefore, it is consistent to replace Eq. (5.29) 
by 



1 + e 


cos 



(5.30) 


If we multiply out the right-hand side of the above equation and 
then express the powers of the cosine in terms of cosines of the 
multiple angles, and neglect term of order r\2 , we have: 



(l + 3e 2 + ^-e 4 ) + e (4 + 3e 2 ) cos + e 2 (3 + cos 2G^ 

- 4 

+ e J cos 3G^ + cos 4G^ (5.31) 


which, on 




3 


integrating, gives: 

2 

(l + 3e 2 + e 4 ) F 2 + e (4 + 3e 2 ) sin G.^ + ^p(6 + e 2 ) sin 2G 1 

e 3 e 4 

+ -j- sin 3G^ + -j 2 " s; *- n 4G q (5.32) 


where, as in Eq. (5.28), we have used the additive constant to 
replace f by F 2 in the secular term. 
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Again noting Eq. (5.7*), we next evaluate the combination 
( cf> 2 “ n 2 4>3 ) . Combining Eqs . (5.28) and (5.32) in this manner 

we have : 



[l + 4n 2 (1 - 2N 2 )] 

en / lc ^ 

*1 

fail ^ / is- j J 


- n 2 e |^[4 + 3e 2 + ^-(7 - 13N 2 )] sin G 1 

+ f [ 6 +e 2 -xf - n2 )] sin 

/ 2 \ 1 - N 2 . 1 

• sin 3G 1 + ^-(e+N -ljsin 4G X + “ 54 — sin 6G x 

(5.33) 


which, using an obvious notation, we may write as 


*2 ~ n *3 = b 0 F 2 +e 


b n sn (nF 1 ,k 1 j - n 2 ^ 3 n sin nG-J . (5.34) 


We now multiply Eq. (5.7*) by the factor 


1 - 


2 „ 2 
n 1 


4 „ 4 

n a 


and then substitute from Eqs. (5.19) and (5.34). After combining 
the secular terms we have: 
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(>- 


2 n 2 4 „ 4 

n & 

2 8 

- n 2 N <|4r tl + 


r-)q 1 ('t’ - a 0 ) = tan* 1 ^ sc(F 2 ,k 2 )] 
[^( X + ^( 3 - h2 )) + s l b o( 3 - 


l( X * 2 -r-) e [S i b n (sn(nF 1 ,k 1 ) - sn(nj lW/ k x )) 

4 nL x (i - nG l ” sln 5 ^; j l“j 

+ ^-(l - N^) sin 2G2 


+ s 


+ n 


(5.35) 


where we have replaced fl o b Y — an adjustment to compensate for 
the constant terms introduced on the right. The angle can now 
be interpreted as the angle of the "first" nodal crossing, that 
is, when f = - 10 , we have c{> = f2 0 . 

We now introduce the final notation for the expression of the 
formula for <p . We set* 



l 

J 3 m p2 32 

so that Eq. (5.35) may be rewritten: 


(5.36a) 
(5.36b) 
(5.36c) 
(5 . 36d) 


*In fact, by making the necessary series expansions we find 
at least to second order that j 3 = 1. 
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j 3 (<l>-«o) = tan_1 [ N sc ( F 2' k 2)] - n2Nj 3 



If we define the angle by the relation 



then the relation [Eq. (5.37)] takes the compact form: 

tan j^j 3 (4> - = N sc(F 2 ,k 2 ) , (5.39) 

which is a clear generalization of the corresponding relation in 
the Kepler problem. 

Note from Eq. (4.11b) that the crossing of the equatorial 
plane (cos a = 0) corresponds to the vanishing of sn(F 2 ,k 2 ) and 
so, from Eq. (5.39), to the vanishing of tan [ j 3 ( - ^) 1 / so 
that 0, represents the angle of the ascending node and the secular 
and periodic variations of £2 can be read off directly from 
Eq. (5.38). 
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6. THE TIME-ANGLE RELATIONSHIP 
To integrate relation (2.26), we first write it in the 

form: 

. dt 1 . 2 2 tc i \ 

A ^ j + c cos o (6.1) 

u 


where we have replaced R by 1/u in accordance with Eq. (3.18) . 
We now introduce u and cos a from Eq. (4.11) and, noting the 
relations (4.10a) and (3.5), we get: 


_A_ dt 

2 df 
P 


1 + 6 cn 1 


1 + e cn| 

( F i' k i). 


+ 


n q Q s 


0 (i + n 2 4a 0 ) 2 (] 




( 6 . 2 ) 


If we integrate Eq. (6.2) we have: 


^i 4( t- fc o) = 

p 


H i + n q Q s. 


o (l+r, 2 e 2 d 0 ) 2 ji( 1 -N 2 )H 2 (6.3) 


where jn is given by relation (4.1) and t Q is the constant intro- 
duced by the integration. H-^ and H 2 are given by 


H 1 = 


f 

1 

+ 

6 

cn 

^ F i ,k i) 


1 

+ 

e 

cn 

( F i' k i)_ 


dF. 


(6.4a) 


H 


2 




dF. 


(6.4b) 


we must evaluate H^ and H 2 individually. 


45 



In considering we first note from straightforward decom- 
position that 


1+6 cn 


( F i- k i) 


1 + e cn^F 1 ,k 1 ) 


(l - n 2 a ) 2 

[l + e cn(F 1 ,k 1 )] 2 

+ . JnW l-. &) . + n 4 a 2 

1 + e cn^F^,k^^ 


(6.5) 


where we use relation (4.12a) for 6. Also it can be readily 
checked that 


jl + e cn(F^,k^jj 2 1 - e 2 ^1 + e cn^F^k^J dn^F^k^) 


d 

e sn 

( F i- k i) 

J 

dF l 

1 + e cn 

( F 1 ' k i) 

J 


( 6 . 6 ) 


so that, combining relations (6.5) and (6.6), we have 

2 


1+6 cn(F^,k^) 


1 + e cn 


( p i' k i) 


1 + e cn^F^kjJ 


, _ 2 , 2 4,2 /. 0 2\ 

1 - 2 r) de - T) d (1 - 2e ) 

1 - e 2 

(l - n 2 d) 2 l 


d 


e 

sn ( F i' k i) 

dF l 

1 

+ 

e cn(F 1 ,k 1 ) 


1 - e' 


, 4,2 

+ n d 


an ( F i ,k i) 


(6.7) 


Accordingly we set: 
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(6.8a) 


L = f dF l 

1 J 1 + e cn(F lf kJ 


L = [ 1 d 

2 J dn(F 1 ,k 1 ) dF l 


e sn 


( F l' k l) 


1 + e cn 


( F i- k i) 


dF. 


(6.8b) 


and, therefore, from relations (6.4a), (6.7), and (6.8) we have: 

H 1 = — 2n 2 de 2 - n 4 d 2 (l - 2e 2 )J 1^ 

- (l - n 2 d) l 2 | + n 4 d 2 F 1 , 


(6.9) 


so that the determination of Hj_ is reduced to the evaluation of 
Lj and L 2 - 

To proceed with the calculation we first observe that if 
we set: 




- e 2 sn(F. ,k. ) 



yi - e 2 

e 

+ cn ( F i' k i) . 


it can be immediately verified that 

dl x dn ( F l' k l) 

dF l 1 + e cn^F^kjJ 


If we also introduce the notation 


I 


12 



sr2 ( F 1 ' k l) dF l 

cn ( F l ' k l )][ 1 + dn(F 1 ,k 1 )J 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 
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then, from relation (6.8a), it follows that 


L, = 


I 


dF, 


1 + e cn 


( F i' k i) 


f dn| 


) dF 1 

+ 

1 . 

r. 

1 - dn| 

l F l' k l)] 

J 1 + e cn ( 

! F i' fc i) 

1 1 + e 

cn | 

( F i- k 


I 1 + k l X 12 


Also we set: 


e sn 


*2 ' 


( F l’ k l) 


dn 


( F i- k 1 ) 


1 + e cn 


( F i ,k i) 


and 


■21 


In 


e sn 2 ^F^ ,kjJ cn ,k-^ dF. 


1 + e cn(F^,k.jJ dn^^F^,k^^ 


Then an integration by parts on relation (6.8b) gives 


L 2 I 2 k l I 21 ' 


If we introduce relations (6.13) and (6.15) into relation ((. 
and neglect terms of order r)6, we have: 


H 1 = 


i-j j 1 - 2n 2 de 2 - n 4 d 2 (l - 2e 2 )] I ± - (l - n 2 d) 2 i. 

+ k 2 [(l - 2n 2 de 2 )l 12 + (l - 2n 2 d)l 21 ] | + n 4 d 2 F 1 


(6.13) 

(6.14a) 

(6.14b) 

(6.15) 

>.9) 
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and so 


1 ifl- 

2n 2 de 2 - n 4 d 2 (1 - 2e 2 )]l 1 - (l - r\ 2 d) I 2 


1 - e 2 ’ 


+ n 2 g 2 e 2 [ 

(l-2n 2 d)(i 12+ i 2l ) + 2n 2 d(l-e 2 )i 12 ]| 


+ n 4 d 2 F x 


(6.16) 


For consistency, we need to calculate (I 12 + I 2l) U P to fi rs t 
order in n 2 , but it will suffice to calculate I 12 to zero-th 
order in r) 2 . 

From relations (6.12) and (6.14b), we have 



| ’sn 2 (r 1 ,K 1 )dF 1 

J dn 2 ( F l' k l) 



e 



1 - e' 


1 + e cn(F 1 ,k 1 )J|dn 2 (F 1 ,k 1 ) l+dn^-^kjj 


dF . 


(6.17) 
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Next we note that 


dn 


( F i' k i) 


= dn 


( F i- k i) 


+ k 


3k 


sn 

dn' 

2 


|( F i- k i) 

( F X' k l) 


1 + 


dn 2 (F 1 ,k 1 ) 

1 + dn(F 1 ,k 1 ^ 


= dn ( F i' k i) + ~T" sn2 ( F i' k i) + °( n ) 

and, alternatively, 

2 ( F i- k i) 


sn 


dn 


( F i ,k i) 


= 1 + k 


1 dn 2 (F 1 ,k 1 ) 

1 + k 2 sn 2 (F 1 ,k 1 ) + o(n 4 ) 


From a similar calculation we note that 


1 + 


dn ( F i' k i) 


dn 


( F l' k l) +k l sn 2 (F 1 ,k 1 )*-t 

v 1 1/ 1 v 1 u |l + dn(F 1 , 


= |[ dn ( F i Fk i) + - r~ sn 2 ( F i' k i) + °( n4 ) 


and, alternatively, 


1 1 

, . .2 sn2 ( F l' k l) 

1 + dn^F-^k-jJ 

+ kl t 1 + d„( F1 , kl jpJ 

II 

nj|h 

1 + -^ sn 2 (F 1 ,k 1 ) + o(n 


(6.18a) 


(6.18b) 


k i)] 

k i)] 2 J 

(6.19a) 


(6.19b) 
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By combining relation (6.18a) with relation (6.19a) and relation 
(6.18b) with relation (6.19b) and omitting terms of order n # we 
get, respectively: 



In the second integral in relation (6.17), we now use relation 
(6.20). With the first term in the first bracket we use rela- 
tion (6.20b), and with the second and third terms we use rela- 
tion (6.20a). If we also use relation (6.18b) in the first 
integral and neglect terms of order we get: 
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If we observe from relation (6.12) that, except for terms of 
order n 2 / we may write: 


12 


1 f sn 2 (Fj, ki)dFi 
( F l' k l) 




(6.22) 


+ e cn 


then it follows that 


d )( I 12 + I 2l) + 21,2,3 

(!-e 2 : 

hia = 

[l- n 2 ( 2 d + 

W)] 

f 


, 2 

F . sn( 

F. 

1 sn 2 (F 1 ,k 1 )dF 1 + (l- 

2n 2 d) 

1 - e 

T 1 + } 

1 11 

2e 2 

I 1 2 + 
x 2e^ 

2e 


+ k 


2 Jsn^(F 1 ,k 1 )aF 1 + n 2 (l-e 2 )(d + | g 2 )j! n(Fl '^ dFl 


Q 2 2 

+ 9r > ?- e 
8 


J sn 2 (F 1 ,k 1 ) cn (F 1 ,k 1 ) dF-j^ 


+ e cn(F 1 ,k 1 ) 

(6.23) 


We next note that except for terms of order n we may write 

r r a rr, c ° s 2g i c ° s 4c ii 

Jsn 4 (F 1 ,k 1 )dF 1 = Jsin G ± dF 1 = J |§ 2 + 8 J dF l 


sin 2G. sin 4G, 
= I F 1 4 + 3T- 


(6.24a) 


J 


sn 2 (F^kJdF, 

1 + e cn (f x /k^ 




1 - e cn(F^,k^J - 


1 - e" 


1 + e cn 


( F i> k i) 


dF. 


F^ sn 
~2 


( F l' k l) 1 - e 2 - 

~~i 2 X 1 


(6.24b) 
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cos dG-^ 

sin G^ sin 3G^ 
4 12 


(6.24c) 


where G^ above is that given by relation (5.25d). If we substi- 
tute relations (6.24) into relation (6.23) and rearrange, we see 
that 


(l-2n 2 a)(l 12 + I 21 ) + 2n 2 a(l-e 2 )l 12 = [l - r, 2 (2d + lg 2 )] 

•j" sn2 ( F l' k l) dF l- ^ij 1 - n 2 [2d(2- e 2 ) + |g 2 (l- e 2 + |e 4 )]J 


+ U - n‘ 


[ 2 a(2-e 2 ) + |g 2 ( 1 -e 2 )]j 


1 - e^ 


sn 


L 2e 


2 1 


In + 


( F l' k l) 


] 


2219 e 3 e 1 

+ n g e yjsin G^ - -^-sin 2G^~ -^sin 3G^+ j^sin 4G.j| (6.25) 


Finally we evaluate the first integral on the right of relation 
(6.25). Noting that 


2 

sn 



we have 


1 - cn(2F 1 ,k 1 ) 

1 + dn(2F-^,kjJ 


(6.26) 



(6.27) 


Using formulas corresponding to relation (6.19b) in the first in- 
tegral and the analog of relation (6.19a) in the second integral. 
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we 


see, after neglecting terms of order n , that 


Jsn 2 (F 1 ,k 1 )dF 1 = | JdF 1 - |Jon(F 1 ,k 1 )an(F 1 ,k 1 )dF 1 

+ -jp [j sin2 2G 1 dG l 

- 3 (sin 2 2G^ cos 2G^ aG^J 


, . k ih F i sn ( 2F i' k i) 

1 T 2 ■ 2 


k iT . 

- 35L 3sln G 1 • 


sin 4G 


sin 3G^ + 


*] 


(6.28) 


If we substitute from relation (6.28) into relation (6.25) , we 
get : 


(l - 2n 2 a)(l 12 + i 21 ) + 2n 2 d (l - e 2 )i 


12 


= 1- n 2 [2d(2 - e 2 ) + |g 2 (l-e 2 ) 


2 \ 1 1 1 - e" 


sn 


2e 


2 X 1 


I, + 


( F i' k J 




. 9 - 7e _2 

+ — q — g 


]] F i 


- [l- n 2 (2d + Jg 2 )Jsn(2F 1 ,k 1 ) 

2 2 r 

+ 2— 2^.1 3 (3 - e) sin G-j^ - 8e sin 2G ± 

- (3-e)sin 3G^ + ^ sin 4G^ , 


(6.29) 
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which we now insert into relation (6.16). To make the resulting 
formula more compact we set: 


D - [de 2 - V(l-e 2 )] + ^[d 2 (l-2e 2 ) + g 2 (l - e 2 ) 

•[d(2-e 2 ) + e 2 )]J , (6.30) 

and insertion of relation (6.29) into relation (6,16) yields: 

Hi = -2n 2 D)i! - (l-n 2 d) 2 i 2 J 

- - n 2 [d 2 ♦ 2g 2 d + ,*]] 

+ ^ [ 2d ( 2 ' * 2 ) + t 9 ^ 1 ' e2 )]] sn ( p l’ k l) 

- ejl - n 2 (2d + Jg 2 )j sn ( 2F i' k i)| 

4 4 3 f 

+ 71 ^ e 2 \ 3 ( 3 - e ) sin - 8e sin 2G^ 

32 (l - e jL 

- (3- e) sin 3G^ + ^ sin 4G^J . (6.31) 

Before we evaluate H 2 we first note that in Eq. (6.3) we 
wish to interpret t Q as the time of first perigee passage. With 
this in mind we have kept free of any constant term, that is, 
satisfies the condition: 


f = 0 =£► H 1 = 0 . 


(6.32) 
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We shall therefore also require that 


f = 0 =*> 



(6.33a) 


or in terms of F 2 , this means 


F 2 - “ =*■ H 2 = 0 


(6.33b) 


In the evaluation of H 2 , we proceed as we did in going from 
relation (6.27) to relation (6.28) and introduce constant terms 
to satisfy relation (6.33). On neglecting terms of order t \^ , we 
get : 

H 2 = f sn2 ( F 2' k 2 ) 3F 2 

§ - |[sn(2F 2 ,k 2 ) - sn(2 M ,k 2 )J 
k 2 

- | 3 Jjsin G 2 - sin y 2 J - £sin 3G 2 - sin 3y 2 J 


= 1 + 


+ j[sin 4G 2 - sin 4y 2 . (6.34) 


Ji- 


where G 2 is as given by relation (5.16c) and y 2 is as given by 


Y 2 



(6.35) 


If we now divide Eq. (6.3) by and introduce Hy and H 2 
from relations (6.31) and (6.34), respectively, we get, after 
combining the secular terms: 
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4(‘-‘o)- 


l fi-2n 2 D I _ (1 - n 2 d) 2 1 

l-e 2 L H 1_ Jl 2 J 
[ vJ 1 '" 2 ) - 9 2 + 2 n 2 [ 


2 - 
n f 


2 £ 2 ( 1 - 
+ e d + - KX 


g 2 + 2n 2 d 2 + 2g 2 d + 9 ^ g 4 




l-n 2 [2d(2-e 2 ) + §g 2 (l-e 2 )] 


sn 


( p i' k i) 


2 2 1 - n 2 (2d + Xg 2 ) 

- n [^2 ' ^ =n( 2F r k !) 

- \^(l + 2n 2 e 2 d 0 )(l-N 2 )[sn(2F 2 ,k 2 )- sn( 2 U ,k 2 )]l 

4 ( 4 3 r 

+ ^ 3 (3-e)sin G, - 8e sin 2G, 

( 1 - e L 

- (3 -e)sin 3G^ + ^ sin 4G J - £ 2 (l - N 2 ) 

• ^3 (sin G^ - sin Y 2) ~ ( s ^- n 2<2 2 - s i n ^ 2 ) 


+ ^sin 4 G 2 - sin 4y, 


(6.36) 


where we have substituted for k .2 from relation (4.12). Before 
attempting to simplify the above, we first consider the multi- 
plying factor on the left. 


We observe that 



(6.37) 
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I 


We use relation (4.10a) to substitute for the first factor on the 
right and relations (3.48) and (3.2) to substitute for the second 
and third factors, respectively, to get 



q s (l 
^o o\ 


+ 


2 2 , 

0 e d 
o o 


r i 


i + 


2 2 

n % m 
o o 



(6.38) 


which, if we substitute for p Q from relation (3.5), reads: 



q s 

^o o 




1 + 


2 2 A \ 2 
0 e d ) 
o O' 


i 


1 + 


2 0 2 

n £ o 



(6.39) 


We note from relation (4.10b) that 


e 


1 + 
1 + 


(6.40) 


and so, after some manipulation. 


1 - e = 


(i - e 2 )(l - n 4 e 2 d 2 ) 

(i * n 2 «&) 2 


and hence 


1 - e' 


(l + "W fl _ 2\ 
= - — -on— l 1 - e ) 


1 -x -I 

' " e o d o 


and so, from relation (4.20), we have: 


t 2 - s (1 - e 2 ) - S ° (l+n ^ } !( 1 - e 2 ) , 
o oV o) i _ 4 2,2 \ / 


1 - 1 e o d o 


(6.41) 


(6.42) 


(6.43) 
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which, with relation (6.39), gives: 


A" 


n . 2 2 , 

1 + n e d 
o o 




If we write: 


o o 


and define a "mean motion" n by setting 


n, V5 


then Eq. (6.44) reads: 


(6.44) 


(6.45) 


(6.46) 


A 


> 2 = (l - e 2 ) 3 / 2 n * 


(6.47) 


If we further define a mean anomaly M by the relation: 


M = n(t - t Q ) , 


(6.48) 


then, by multiplying Eq. (6.36) by the factor 


(i - £ 


213/2 


we obtain the following relation between the two anomalies M and 
f, namely: 
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This formula can be made a little simpler by using the algebraic 
formulas for the constants as developed in sections 3 and 4 and 
by exploiting the fact that we neglect terms of order We 

also set: 
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. _ l - 2n d 

34 ' 3i3t ’ 


3q = 


(l - n 2 d) 2 


3i3 


1 J T 


(6.49) 


and write Ip and I 2 explicitly. After some manipulation with the 
above coefficients, we finally obtain: 


M = j arctan 


+ n 2 e 


f - e 2 sn(j 1 f,k 1 ) 

L e + cn(j 1 f,k 1 ) J 

3/2 

1/2 


-3, 


1 

e sn 

(h f ’ k i) 

V' k l) 

1 + e 

cn ( j l f ' k l) 

)n 2 (21 - 

37N 2 ) 

]‘ 


• 3N 2 + 69N 4 


(i-- 2 ) /2 [(i-n 2 ) - ^[ 41-9= 

- e 2 (l + 58N 2 -75N 4 )]Jsn(j 1 f,k 1 ) + n 2 (l-e 2 ) 1/2 

• | e2 ^(!-N 2 ) - ^[39 - 122N 2 + 99N 4 + e 2 (8 - 76N 2 + 68N 4 )JJ 

• sn(2j 1 f,k 1 ) - H 1 ~ N ^ l + n 2 |^5N 2 + 3e 2 (3N 2 - l)]j 

• |^sn [2 (f + u>)k 2 ] - sn(2o),k 2 )j|+ ^(l-N 2 )(l- e 2 ) V2 

• |(l - N 2 ) e 3 £3 (3 - e) sin - 8e sin 2G 1 - (3 - e) sin 3G 1 

2 

+ | sin 4G 1 J - (l-e 2 ) £3 (sin G 2 - sin y 2 ) 

- (sin 3G 2 - sin 3y 2 ) + | (sin 4G 2 - sin 4y 2 )]j . 


(6.50) 


When n = 0, only the first two terms survive in which case the 
above relation reduces to the well known time-angle relation of 
the Kepler problem. 
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7 . APPROXIMATE FORMULAS 

Equations (4.11) and (5.39), with ft given by Eq. (5.38), 
give the formulas for the three coordinates in terms of f. These 
are complemented by the time-angle relationship (6.50). 

In arriving at relations (5.38), (5.39) and (6.50), we have 

made second-order approximations. It is therefore consistent to 
introduce such approximations into all the relations. Though this 
means that the formulas will, to a certain extent, lose their 
compact form, they will become simpler in that they will no 
longer involve elliptic functions. The latter will be replaced 
by their second-order approximations in terms of truncated 
trigonometric series. For clarification we list the sequence 
of steps leading to these approximate formulas; the analysis and 
derivation are available in the standard treatises, e.g., Whit- 
taker and Watson (ref. 9) or Davis (ref. 10) . 

For the modulus k of the elliptic functions appearing in 
either Eqs. (4.11a) or (4.11b), we note that 

k 2 = o(n 2 ) , (7.1) 


and, if we define the complementary modulus k' by 


k' 




/ 


(7.2) 


then we have that 


^=( i -* 2 ) V4 = 1 -£-^ 4 + 0 („ 6 ) 

1 - VP = !j-|l + §k 2 + ^k 4 + o(n 6 )j (7.3) 

- ^k 4 + o(n 6 )] 
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Defining e by the relation 


- = 1 1 ~ y/EL 
2 1 + 


(7.4) 


then, after expanding, we find that 


c k 2 f, A k 2 x 21,4 X ./ 6\1 _/ 2\ 

e = 16 L 1 + I" + 64 k + °( n )J = °( r ' ) * 


When the modulus q of the associated Theta functions is written 
in terms of e , we have 


q » e[l + 0 ( e 4 ) J 

- £[l ♦ 4 + ^ + °(^)] • ( 7 - 6 ) 

and so we have: 

1/2 k I", , k 2 17 .4 ~( 6 \"1 

g = 4L 1 + t + i 28 k + °( n jj 

( 7 . 7 ) 

qi/2 _ x T k 2 17 . 4 / 6\~| 

k 4 L 1 + 4 + 128 k + °( rt )J 

The period 4K is related to 2tt as follows: 

(f) V2 = 1 + 2 q + o(r, 6 ) 
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and so we have 


2K 

TT 





TT 

2K 


= 1 




+ o(n 6 ) 


(7.9) 


If the angle G is then defined in terms of F by the relation 


G = 






t 


(7.10) 


then the trigonometric series representations of the Jacobian 
elliptic functions take the form: 


„ * ^ n+1/2 

sn (f ,k) = kJ 2-/T 3 : 2n+l sin (2n + 1)G 

n=0 ± q 


1/2 

= 4 a - 11 


E 


n 


k 2K i n 2n+l 
A 1 - q 
n=0 


sin (2n + 1)G 


oo 


n 


n=0 1 - q 


2n+l 


sin (2n+ 1)G (7.11a) 


cn 


4 00 

(F,k) = (l ~ Y28 + **•) S 


n 


A 1 - q 
n=0 


2n+l 


cos (2n+ 1)G (7.11b) 


dn ( F , k ) = (l- T + ") 1+ 4 1+ ^2n+l 

L n=0 q 


] 


cos 2nG I (7 . 11c) 
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oo 


am 


(F,k) = G + ]T 


2q 


n 


( 1 2n\ 

n=l "l 1 + q ) 


cos 2nG . 


(7 


We also include the trigonometric series representation of 
sc (F ,k) , namely : 


>c(f ,k) = ^ 


! + id + ii *1 + 

4 64 


[ 


tan G + 4 


E f v n 2n 

sin 2nG 

n=l ^ 


]■ 


(7 


In a second-order approximation the above series can be 
truncated. The approximate forms we shall use are: 


sn(F,k) = 


( 1 + t? + It?) 

* £(i - 4 ) 


sin G 


sin 3G + 


256 


sin 5G 


(7 


en(F,k) = 


, k 2 9k 4 
16 256 


cos G 


k 2 (. k 2 \ . k 4 

+ 16 V 1 + 2 ) cos 3G + 256 


cos 5G 


(7 


dn(F,k) = 


. * 2 . 

1 

M 

LO 

+ k2 (i 

/ T 
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-> 


1 + -|-)cos 2G 


+ cos 4G 


(7 


SC 


- (i + 4 + #) 


tan G - g-j sin 2G 


(7 


.lid) 

.lie) 

• 12a) 

• 12b) 

.12c) 
. 12d) 
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with 


G 




(7 . 12e) 


We shall now introduce these approximations into the formulas for 
the three coordinates . 

We first note that, in accordance with relation (7.12e), we 

have : 


G 


1 



1 - 


k 2 

k l 

"T 



(7.13a) 


G 


2 





to) 


(7.13b) 


These relations have already been anticipated in sections 5 and 6 . 
We next introduce the above approximations into Eq. (4.11). If in 
Eq. (4.11a) we also substitute for 6 from relation (4.12c), we 
obtain 



(7.14a) 


and 


cos a 


- N 2 





sin G 2 


+ 




sin 


3G. 



sin 



(7.14b) 
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If we also introduce the approximation into Eq. (5. 38' for ft, we 
have : 


ft = 


ft Q - r| 2 ni g N (f + co ) - r) 2 fl + ^y m pl Ne 


r / \ ( . . ^i* \“| 

^sin - sin cjj + b 2 ^sin 2G 2 - sin — u)H 

jl e [ ZX( sin nG l ' sln _ 2iq i “)] 


+ n 4 N 


+ m 


P 2 


(l - N 2 ) 


sin 2G, 


in which the coefficients B n are given by 


B 3 = 3 3 - 


e 2 (i - N 2 )b x 


16 


' B 6 ~ 3 6 


( x - h2 )' 
16 


(7.15) 


(7.16) 


B =8 for n = 1, 2, 4, 5 
n n 


The approximate formula for the third coordinate is obtained by 
introducing relation (7.12d) into Eq. (5.39). We get: 



With ft as given by Eq. (7.15) and relations (7.16), then 
relations (7.14) and (7.17) are the second-order approximate 
formulae for the three coordinates in terms of the "true anomaly" 
f. We can get a corresponding form for the time-angle relation- 
ship by introducing the approximate forms into Eq. (6.50). Be- 
cause of the unwieldy length of the resulting formula, we do not 
exhibit it at this time. 
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8. A GEOMETRICAL RESULT 


In the case of the Kepler problem, the conservation of 
angular momentum implies that the motion of the particle is 
planar, and it can be immediately shown that for negative energy 
the particle moves between two concentric circles in the plane 
of motion. In this section, we seek the analog for the Vinti 
problem of the above geometrical result. 

Here the motion is no longer planar and, in general, is 
quite complex. However, if we assume that the exact formulas 
(4.11) and the approximate formulas (5.38) and (5.39) for the 
three coordinates indeed give an exact discription of the motion, 
we can get some geometrical insight. 

To describe this, it is again useful to think of the Vinti 
problem as a perturbation of the Kepler problem. If we give the 
plane of the corresponding Kepler problem an appropriate deforma- 
tion and then allow this surface to rotate with the angular 
velocity of the instantaneous nodal line (S') derived from Eq. 
(5.38), we can show that the particle remains on this surface. 

More precisely, we shall show that the motion takes place 
in a torsidal region defined by two ellipsoids of revolution and 
a hyperboloid of revolution. In this region the particle remains 
on a surface S rotating with the angular velocity of the instan- 
taneous node, the surface S being a "small" deformation of a 
plane. The remainder of this section is devoted mainly to the 
explicit derivation of this result. 

We first derive a relation between a and (j) . If we set: 


^ = j 3 (4> - 8) , (8.1) 

then Eq. (5.39) takes the form: 

tan t|> = N sc(F 2 ,k 2 ^ . (8.2) 

It follows that 

tan^ ip = sc^^F 2 ,k 2 ^ , (8.3) 
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and so 


sec 


2 


Also, from Eq. 

sin 2 


Combining Eqs. 


2 

sec 


4> = 1 + N 2 sc 2 ^F2,k2^ 

cn 2 (F 2 ,k 2 ) + N sn 2 (F 2 ,k 2 ) 

cn2 ^ F 2,k 2 ) 

(4.11b), it follows that 
a - 1 - (l-N 2 ) sn 2 (F 2 ,k 2 ) 

= cn 2 ^ F 2/k2^ + N 2 sn 2 ^ F 2,k2^ . 

(8.4) and (8.5), we see that 



(8.4) 


(8.5) 


( 8 . 6 ) 


and hence 


cos j ^ (<j> - ft) = cos ip = 


cn 


(f 2 J 


(8.7) 


sm a 


We introduce a system of rotating Cartesian axes by the 
relations : 


X = 

(r 2 + 

c ) sm a cos 

(<P~ 

ft) 



Y = 

(r 2 + 

2\^/^ 

c ) sm a sm 

(<t> - 

ft) 

1 

f (8.8) 


Z = R cos a J 
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Ill 


i mi I 


II nil 


IBII III 


where we note that X = 0 corresponds to <|> = fi. In relation to 
the fixed x-y-z system of Eq. (1.1), the Z-axis of Eq. ( 8 . 8 ) 
coincides with the z-axis of Eq. (1.1), while the X-Y system 
rotates at the angular velocity [obtainable from Eq. (5.38)] 
relative to the x-y system. Thus the X-axis always coincides 
with the instantaneous nodal line. 

We further introduce an auxiliary system of variables (£]_, 
5 2 / £ 3 ) defined by 


= R sin a cos 
£ 2 = R sin a sin 
£3 = R cos a 


j 3 (<j> - ft) = R sin a 
j 3 (<|) - ft) = R sin a 


cos \p 

sin \p (8,9) 


and it immediately follows that 


^2 

? 1 


tan ip = N sc 




( 8 . 10 ) 


where we have used Eq. (8.2). From Eq. (8.10), we have 
h = H . Vl - N 2 sn(F 2 ,k 2 ) 

? 1 V 1 n2 cn ( p 2' k 2) 


which, if we use Eq. (4.11b) in the numerator and Eq. (8.7) in 
the denominator, gives: 


^2 

h 


N 


cos a 


_ n 2 sin a cos \p 

?• 

R 


N 


yi - N 2 


( 8 . 11 ) 
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the latter following directly from Eq. (8.9). If we define an 
"asymptotic” angle of inclination I by the relation 


cos I = N , 


( 8 . 12 ) 


then Eq. (8.11) takes the form: 


- C2* tan 1 • 


(8.13) 


Using the above with Eqs. (8.8) and (8.9), we have that 


Z = = ? 2 ' tan 1 


= R»sin a*sin j ^ ( ^ -fi)*tan I 



+ *c 2 ) 1/2 



+ 



sin a*sin 


(<f> - «) 


sin j _ (<p - ft) 

- • tan I 

sin (<}> - fi) 



R 


+ 



sin j ^ (<j) - ft) 

•Y - tan I . 

sin (<J> - ft) 


(8.14) 


If we set j-j = 1 + e , so that e = 0(n 2 )* and use the addition 
formula for the sine term in the numerator of Eq. (8.14), we get: 


Z = Y • tan I- 


( p 2 , _ 2\ 1/2 


R + cj 

' [cos e (<j> - ft) + cot ((J> - ft) sin"e((|>-ft)] . (8.15) 


*As remarked in the footnote on page 43, we have, to second 
order, j 3 =l,i.e.,r=0(n 6 ). 
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Again using Eg. (8.8) to substitute for cot (<f> - R) , we have: 


Z = tan 


( r2 + c 2 ) 


1/2 


Y cos e (cf> - R) + X sin e (cf> - R) . (8.16) 


When r| = 0 , we also have R' = 0 [see Eq. (5.38)] and the 
above equation takes the form 


Z = Y • tan I , 


(8.16*) 


which is the equation of the plane of the Kepler motion. For 
ri 0, Eq. (8.16) represents a surface, clearly a deformed plane, 
moving with angular velocity R' relative to fixed x-y-z axes. 

We let S denote this surface. 

As a description of S, Eq. (8.16) is incomplete since we 
have not expressed the quantities R and (<j> - R) in terms of X, Y 
and Z. This involves the solution of the following set of trans- 
cendental equations, obtained directly from Eq. (8.8), namely: 


R 2 + c 2 sin 2 a = X 2 + Y 2 + Z 2 (8.17a) 

R cos a = Z (8.17b) 


and 


tan (cf) - R) = Y/X . (8.18) 


An adequate approximate solution of the pair of Eqs. (8.17) 
can be obtained by expanding in powers of (c/R) 2 to second order. 
However, since there is little to be gained, we do not include 
this calculation here. We only note that Eqs. (8.17) and (8.18) 
together complement Eq. (8.16) to give a complete description of 
the surface S. 

It is clear from relation (3.17) of section 3 that the motion 
takes place between the two cenfocal ellipsoids of revolution 
defined by R = R^ and R = R 2 where R^ and R 2 are the roots of the 
first quadratic on the right of Eq. (3.17). We note furthur from 
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Eq. (4.11b) 


and using Eq. (8.12) 


that 


cos a 


L Vl - 


N = sin I 


(8.19) 


so that the trajectory is bounded by the hyperboloid of revolu- 
tion : 


cos a | = sin I 


( 8 . 20 ) 


Hence the motion takes place in a toroidal region consisting of 
that portion of the "exterior" of the hyperboloid [Eq. (8.20)] 
lying between the ellipsoids of revolution defined by R = and 

R = R 2 . 

We can now give a geometrical interpretation to what we 
termed the asymptotic angle of inclination I. The bounding 
hyperboloid [Eq. (8.20)] has an asymptotic cone: the angle I is 

the angle of inclination of the tangent plane to this asymptotic 
cone . 


9 . THE POTENTIAL 

By an elegant device involving the introduction of complex 
variables and the use of the generating function for Legendre 
polynomials, Vinti (ref. 1) has shown that when the potential 
[Eq. (1.5)] is expressed in terms of spherical coordinates r, 9, 
and <p , we get: 


l / 


y R 


R +c 


cos 


= -£ I>>"(?) 2n ' 


P 2n (c°s 


0 ) 


(9.1) 


n=0 


where the P 2 n term denotes the Legendre polynomials of even order. 
If r Q denotes a length scale, which we shall later identify with 
the radius of the Earth, we can write Eq. (9.1) as 


, x 2n/r \^ n 

-£2>> n (a { 4 } ■ 


P 2n (cos 


6) 


n=0 


(9.2) 
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If we now make the identification 


J 


2 



(9.3) 


then Eq. (9.3) becomes: 


1 / 


y 

r 


°° 2n 

p 2n 


n=0 


(cos 


0 ) 



- E ( -> n+1 

n=l 




(9.4a) 


(9.4b) 


The standard representation of the geopotential in spherical 
coordinates — referred to an origin at the center of mass--is 


( — n 

■'e ■ -r I 1 - E J n(rO P n (c ° s e) 

( n=2 

oo oo 21 

- E E J nm(r) P n (cos 6) cos 

n=2 m=l 


In Eq. (9.5) r Q denotes the mean Earth radius and the J ' s are 
called the geopotential coefficients. 

Apart from the Kepler term the dominant term in the repre- 
sentation [Eq. (9.5)] is that with coefficient J 2 , which exceeds 
by an order of magnitude the effect of any other term. If, for 
purposes of comparison, we restrict our attention to that part of 
the geopotential which is both rotationally symmetric and sym- 
metric about the equatorial plane, then we include this dominant 
term. We denote this part of the geopotential by t/gg and the 
residual (l/ G - I^qs) can b e considered a perturbation which must 
be taken into account later. Then: 


(<P ~ a mn ) 


(9.5) 


74 



u 

r 


(9.6) 


GS 


= n 


oo 2n 

-E J 2 n(^) 1 


P 2n^ COS 


m=l 


We now see that when c is chosen in accordance with the identi- 
fication [Eq. (9.3)], then the Vinti potential [Eq. (9.4)] agrees 
with the symmetric geopotential [Eq. (9.6)] up to the second 
(zonal) harmonic. Further identification would require that 


J 2n = (-) n+1 J 2 for n > 2. 


This is not true for any n > 1 and is not even true in an order 
of magnitude sense for any n > 2. In fact: 

J 4 ~ “I J 2 ' l J 2 J >5> i J 2 1 n for n * 3 ' (9 * 7) 


so that if the Vinti problem is taken as the base solution, we 
must add as a perturbation the residual (1 /qs “ which starts 
with a term in the fourth harmonic with coefficient of order J 2 . 

As far as we know, there has not been proposed a real physi- 
cal situation giving rise to the Vinti potential field. Although 
this question has no relevance for the dynamical problem and its 
relation to satellite orbit prediction, nevertheless it is of 
interest when we consider how closely the Vinti potential approxi- 
mates the geopotential. The latter feature suggests that there 
is a hypothetical geoid whose potential matches the geopotential 
exactly. This is indeed the case. We shall show that a solid 
sphere with an appropriate interior mass distribution induces an 
external potential of the Vinti type. 

We could start by posing the larger question, namely, what 
mass distribution in the geoid gives rise to the geopotential 1 /q? 
Although the restriction is not necessary, we shall here restrict 
our consideration of this question to the simpler case of the 
symmetric geopotential I/qs* Our problem then is to determine the 
mass distribution inside a sphere consistent with an arbitrary 
external field of the form (/qs • 

Referred to a spherical coordinate system we let P(r,0,cj)) 
represent an arbitrary exterior point and Q (r ' , 0 ' , <j> ' ) an arbitrary 
interior point at which is situated the mass element dm. We 
denote the position vector of P and Q by r and r', respectively. 
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and let x denote the angle between r and r'. If we let T(r',0', 

<p ' ) denote the mass density at any interior point, we see that 

dm = x-r' 2 sin 6' dr ' *d0 ' -d<J> ' . (9.7*) 


If we let y denote the gravitational constant, then the 
potential at P due to the mass element dm at Q is given by 



dm 



(9.8a) 



} + Ot) 2 - 2 (r) cos *] 


1/2 


(9.8b) 


where, in deriving Eq. (9.8b), we have used the cosine law. 
Noting that on the right of Eq. (9.8b) we have the generating 
function for Legendre polynomials, it follows that 


— T _ dm 

au Q = -y — 


12(ir) p n (oos X) 


n=0 


= - > XI -1ST 


r ' n P (cos x) 
n 


(9.9) 


n=0 


We now use the addition theorem for Legendre functions (ref. 9), 
namely, when 


cos x = cos 0 cos 0' + sin 0 sin 0' cos (<J>-<f>') (9.10) 


(as is the case here) , we have the following relation for the 
Legendre polynomials : 
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P n (cos x) 


P (cos 0 ) P (cos 0 ' ) 
n n 


n 


+ 2 y] - |^ + ^y -r P™ (cos 0)P™ (cos 0 ' ) cos m (<f> - <j> ' ) . 

(9.11) 


m=l 


Introducing Eq. (9.11) into Eq. (9.9) and substituting for dm 
from Eq. (9.7), we get: 


dU Q = ^r[ r ' np n (cos 0)P n (coS e ')] 

( n=0 r 

oo 

’ir' 2 sin 0 1 dr 1 d0 ' dcf) 1 + 2 ^ n ~^ -j 

z ; r n 
n=l 

r n 

• r 1 n (n + mp p ^ cos 0 ) p ^ cos 0 ' ) cos ' ) 

m=l 

2 I 

•rr' sin 0' dr' d0 ' dcj> ' V . (9.12) 


To determine the potential at P due to the sphere we must 
perform the integration 


u = /du Q 


when the integration is taken over the sphere, namely, r' ranges 
from 0 to r Q , 0' from 0 to it, and <p ’ from 0 to 2 tt. If we take x 
to be rotationally symmetric, i.e., independent of <p ' , then the 
4) ' -integration in the first summation integrates to 2tt, while the 
* -dependence in the second summation integrates to zero. After 
performing the (^-integration, we therefore have: 


77 



I I 


^ p (cos e) f r o r 71 

“ - - 2 - L “^n+1 J I 

n=0 *'0 • / 0 


. n+2 


r'"" J P (cos 0') sin 0'*x*dr' d0 ' . 
n 


0 •'o 


(9.13) 


In all the above manipulations we have interchanged the integra- 
tion and summation operations . 

If we further assume that x is not dependent on r', then the 
r ' -integration in Eq. (9.13) is immediate and we get 


U = 


9 ^/ r o\ n+1 P n (cOS 0) f 77 

- 2 ^ r 0 S(v) (n + 3) J T * P n ( 

n=0 *^0 


cos 0 ' ) sin 0 ' d0 ' . 


(9.14) 


If we now represent the dependence of x on 0 ' in the form: 


x = x < 1 + 


22 t sl p i (cos 0,) 


(9.15) 


£=0 


so that x D denotes the (constant) mean density and the xj^'s 
(£ ^-1) are the dimensionless coefficients for the higher moments, 
then using the orthogonality of the Legendre polynomials we have: 


J n 


(cos 0 1 ) sin 0'*x*d0' = 




X n I P n (COS 0') 

r o 


[■*£- 

L 1 = 1 


)j si 


x^P^(cos 0 ' ) | sin 0' d0 ' 


2x for n=0 
o 




(9.16) 


for n >. 1 
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and substituting Eq. (9.16) into Eq. (9.14), we get: 


n+1 


U = 


- 4 ^ r oMlT I + S (n + 3) Un+ -rr(^) P n (COS 0) ! 


n=l 


4irr 2 1 3 x /r \ n ) 

- -is 3 "' o i + Ew ^hLu It) p „ (c ° s 9 > • (9 - 17) 


n=l 


If we now call the sphere with constant density t d the mean 
sphere and denote its mass by y Q , then clearly 


4irr~ 


o 

3 T o ' 


(9.18) 


and setting y = y^o (the normalized gravitational constant) , we 
see that Eq. (9.17) may be written: 


U 



1 + 


S (n 
n=l 


3t 


n 


+ 3) (2n + 


lY^rf 1 


P n (cos 



(9.19) 


It is clear from Eq. (9.19) that the x n can now be chosen to fit 
an arbitrary axisymmetric potential. We can make U symmetric 
about the equatorial plane by requiring all odd coefficients to 
vanish, that is, by setting: 


T., ,, = 0 for all k , 
2k+l 


(9.20) 


We then get the symmetric part of U which we call U s in the form: 


( °° 3t r 2n ) 

U s = "r I 1 + S ( 2n + 3 )T4n"+ T)(~T ) P 2n (cos 0) [ * (9 ‘ 21) 


n=l 
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Then, comparing Eq. (9.6) with Eq. (9.21), we see that U s be- 
comes identical with l/gg if we take : 


_ _ (2n + 3) ( 4n + 1) _ 

x 2n ” 3 u 2n 


(9.22) 


and all odd coefficients zero. In particular, we can identify 
U s with the Vinti potential 1/ given by Eq . (9.4b) if again we 

have the odd coefficients in the density distribution identically 
zero and the even ones as given by 


_ _ (2n + 3) ( 4n + 1) , , ,n 

T 2n 3 ( J 2 } 


(9.23) 


We have thus produced a physical realization of the Vinti poten- 
tial field. 

We have also produced a density distribution whose potential 
can be matched to the symmetric part of the geopotential. This 
suggests a procedure which, by successive refinement, may lead to 
an approximation for the density distribution of the Earth, which, 
except for layers of sharp discontinuity, should give some in- 
sight into the actual distribution. The above analysis assumes 
continuity — and, in fact, analyticity — in the variables for the 
density distribution. 


10. CONCLUSION 

We have derived a solution of the Vinti dynamical problem in 
the relatively compact form given by relations (4.11), (5.38), 

and (5.39). These relations are clear generalizations of the 
solution of the Kepler problem in terms of true anomaly. As an 
alternative we have the more elementary expanded form in relations 
(7.14) to (7.17). Each form is complemented by the time-angle 
relationship (6.50). 

We have also derived some qualitative results on the motion 
(section 8) and a physical interpretation of the Vinti potential 
in terms of mass distribution (section 9) . 
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^G 

74 

G 

64 

L 1' L 2 

48 

^GS 

75 

G i 

39 

m 

o 

25 

V 

7 

G 2 

35 

m f f m n 
s' Pi P2 

43 

V 1' V 2 

8 

h 

15 

M 

2,59 

w 

19 

o 

h l 

20 

n 

2,59 

x-y-z 

4 

h 

22 

N 

29 

X-Y-Z 

69 


y 14 


82 



Greek 


Symbol 

Page 

Symbol 

Page 

Symbol 

Page 

a 

ii 

e ' 

75 

T 

16 ,13 


42 

X 1' X 2 

10,11 


1,4 


56 

X 3 

6 


75 

Y 

76 

X 

11 


32 

r 

3 

A 

12 

X 

76 

6 

21 

% 

79 

1j> 

68 

e 

63 

u 

2,5 

V^2 

3 

£ 

71 

V 

14 

w l' w 2 

28 

c 

26 


8 

00 

28 

n 

14 

^1~^2~^3 

70 

n 

o 

43 

0 

1,4 

a 

4 

ft * 

o 

32 





ft 

44 


(•) = differentiation w.r.t. time (t) , p. 6 
(') = differentiation w.r.t. anomaly (f ) , p. 12 
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